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Abstract 

We study a generalized 8-vertex model where the vertices are cou- 
pled to a locally varying field. We rewrite the partition function as 
an integral over Grassmann variables. In this form it is possible to 
explicitly evaluate all terms of the hopping expansion. Applications 
of the resulting formula, in particular its relation to 2-D lattice field 
theories with fermions are discussed. 
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1 Introduction 



Rewriting physical systems in different representations is a powerful tool in 
theoretical physics. Alternative representations of a system can stress dif- 
ferent aspects and allow new insights into its physical behavior. Sometimes 
a model is more accessible to numerical methods in a rewritten form or can 
even be solved explicitly. 

Here we study a generalized version of the 8- vertex model ]l|, ||, ||| , where 
now the vertices couple to an external field. This generalized model encom- 
passes a large class of physically interesting systems such as spin models with 
locally varying coupling, polymers in an external field and in particular 2-D 
lattice field theories with fermions. The model first will be represented as 
an integral over Grassmann variables and in this representation one then 
can explicitly evaluate all terms of the hopping expansion. The resulting 
expression is a new representation for the partition function of the original 
model in terms of loops. In this new form it is possible to make connection 
to 2-dimensional lattice field theories with Wilson fermions Q . The formula 
e.g. can be used to considerably simplify the hopping expansion of the lat- 
tice fermion determinant in an external scalar field, since it establishes that 
a large class of contributions (all terms with multiply occupied links) gets 
cancelled. Furthermore, our formula allows to explicitly integrate out the 
external field and to obtain simple loop representations for e.g. the Gross- 
Neveu model (see Section 5 for a brief discussion of these applications) . 

The article is organized as follows. In the next section we introduce the 
model and discuss its use for polymer systems and 2-D statistical models. 
Section 3 contains the reformulation of the partition function as an integral 
over Grassmann variables and the hopping expansion. In Section 4 the traces 
over the hopping generators are computed and the final expression for the 
hopping expansion is presented. The article closes with a brief discussion of 
the properties and applications of our expansion formula. 



2 Definition of the generalized 8-vertex model 

We analyze a generalized version of the 8-vertex model where now the ver- 
tices are coupled to a locally varying external field B^(x). The standard 
8-vertex model @, g § can be viewed as a model of 8 quadratic tiles (ver- 
tices) and each of them is assigned a weight Wi (i = 1, ...8) (compare Fig. Q). 
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Figure 1: The eight vertices (tiles) and their weights Wi. 

Consider now a lattice A which we assume to be a finite, rectangular piece 
of 2Z 2 (the generalization to e.g. a torus is straightforward). A tiling of this 
lattice is a covering of A with the tiles such that on each site of A we place 
one of our tiles with the centers of the tiles sitting on the sites. The set T 
of admissible tilings is given by those arrangements of tiles where the black 
lines on the tiles never have an open end (for sites at the boundary this 
implies that not all 8 tiles can be used there). The partition function of the 
standard 8-vertex model is the sum over all admissible tilings t G T and 
the Boltzmann weight for a particular tiling t is given by the product of the 
weights Wi for all tiles used in the tiling t. 

In our generalization of the model we now add an additional, local struc- 
ture. To each link (x, fi), fx = 1, 2 of A we assign a real- or complex- valued 
field Bfj,(x), ijl = 1,2. The contour C(t) of a tiling t is defined to be the 
set of all links of A which are occupied by black lines from the tiles. Since 
we allow only admissible tilings, each link is either occupied or empty and 
'half-occupied' links do not occur. The partition function of the generalized 
model is now given by 

Z = £ UvH 1 *® II M*)- (!) 

ter i=i ( x ,n)eC(t) 

By ni{t) we denote the abundance of tile Nr. i in a given tiling t. Note 
that the contour C(t) is simply the set of links occupied by black lines and 
does not have an orientation. So in case the field is chosen complex, a 
link variable B^(x) is always counted as it is, and no complex conjugation 
is implied. 
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A second (less general) form of the model is obtained by assigning a 
scalar field tp(x) to all sites of A and setting 

J3 M (s) = yVO)\/W + A) , 

where /i denotes the unit vector in direction [i. In case one of the 99 is 
negative becomes complex, a case which we explicitly included above. 
The partition function now reads 

z = e n^ ni(t) n ( 2 ) 

ter «=i xes(t) 

Here 5(i) denotes the set of all sites occupied by the tiling t. When a site x 
is occupied by tile Nr. 2, this site is counted twice giving a factor tp(x) 2 . In 
case x is occupied by tile Nr. 1, x ^ 5(4) and the factor is 1. For all other 
tiles x is counted once and the factor is f(x). 

The generalized 8- vertex model in its two forms ((!]) , (j2|) encompasses sev- 
eral interesting physical systems (for a detailed discussion of the correspond- 
ing models at trivial external field see || and references therein). For exam- 
ple by setting the parameters in formulation (||) to w\ = W2 = ... = ^8 = 1 
and the external fields to B^{x) = e - 2 J^( x )/ kT ; the model describes the Ising 
model with locally varying coupling Ju(x) for links (x,fx). The set C(t) has 
the interpretation of a Peierls contour which separates patches of up and 
down spins. Ising type models with next to nearest neighbor terms can be 
obtained by chosing different values for the weights itfj. The model with 
locally varying couplings can furthermore be viewed as a generating func- 
tional for 2n-point functions in Ising type models: Differentiating Z in the 
form (H|) with respect to the link variables on a connected path produces 
the 2-point function for two spins sitting at the endpoints of the path (for 
2n-point functions use a net of paths). 

In its form (||), the generalized vertex model can be used to describe the 
physics of loop gases and polymers (compare e.g. [||, ||, 0, § ||, |l0|) in an 
external field (p. 

Most remarkable, however, is the fact that our model is related to 2-D 
lattice field theories with Wilson fermions. For staggered fermions it has long 
been known, that a polymer representation can be found 
Due to their spinor structure, a similar map for Wilson fermions ]15[ is 
more complicated. Based on the hopping expansion which will be discussed 
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in the next section, it was shown in [Q] that 2-D Wilson fermions in a scalar 
background field ip(x) are equivalent to our model in its form (||) (see Section 
5 for a brief discussion of the argument used in || ) . 



3 Grassmann representation and hopping expan- 
sion 

The first step in computing the hopping expansion is finding a represen- 
tation of the partition function through a Grassmann integral. For the 



standard 8-vertex model this problem has been solved [16, 17, ffg]. The 
central idea is to assign to each site x of the lattice 4 Grassmann variables 
r] + i(x) , r]^i(x) , r] + 2(x) , 77-2(3;) and integrate a Boltzmann factor with an ap- 
propriately chosen action S[rj] to obtain 

Z = (-1)I A I j Y[ dr]- 2 (x)dri + 2(x) dr ] ^ 1 {x)dr ]+1 {x)e s ^. (3) 

Here |A| denotes the size of the lattice. The action is a quadratic form in the 
rji and contains two types of terms. Hopping terms which are products of 
Grassmann variables at neighboring sites create link elements, i.e. the black 
lines on our tiles. In addition, quadratic terms based on only one site are 
needed to saturate the Grassmann integral. These terms enable the different 
moves the lines perform on the different tiles, e.g. a corner on tiles 5-8, or 
keep going straight on tiles 3,4. For an explicit discussion of the Grassmann 
representation for the standard 8-vertex model see fTfJ. 

For the generalized model in the external field B^(x) only the hopping 
terms have to be modified such that the links are furnished with their cor- 
responding fields. The action S = Sh + S m + S c producing the generalized 
8-vertex model then consists of hopping, monomer and corner terms 

Sh = VW2YI [Bi(x) rj +1 (x)r]_i(x + 1) 

x-eA 

+ B 2 (x) T] +2 (x)r]- 2 (x + 2) , 
Sm = —Y]\w3T]-i(x)T] + i(x)+W4 : r)-2(x)T] + 2(x)\, 

[105 rj +1 (x)r]_ 2 (x) + w e rj + 2{x)rj_i{x) 



™2 , eA 



+ W7 r)- 2 (x)r]-i(x) + w 8 r? +2 (x)?? + i(x) . (4) 
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The boundary conditions are open, i.e. hopping terms that would lead out 
of our rectangular lattice are omitted. Inserting (Q) into the path integral, 
expanding the exponent termwise and integrating the fields along the lines 
of |l6j establishes the representation (|3|) for the partition function ([!]). 

We remark, that the action (|j) is independent of w±. However, it turns 
out, that for the Grassmann representation @, (|j) with a quadratic action 
to work, the weights Wi have to obey the free fermion condition Q 

W\W2 + VJ3VJ4 = W5WQ + W7W8 . (5) 

Thus w\ can be computed from the independent parameters Wi,i > 1. Note 
that for the case of a trivial external field (all B^(x) = 1), the free fermion 
condition is a sufficient condition for finding an explicit solution. This is 
most easily seen by first constructing the Grassmann representation (which 
always can be done when (||) holds) and then performing a Fourier transfor- 
mation of the Grassmann variables. The Fourier transformation diagonalizes 
the action and the free energy can be readily computed. However, when the 
fields Bu(x) are chosen non-trivially, Fourier transformation produces a con- 
volution of the external field and does no longer diagonalize the action. 



Here we analyze the hopping expansion for the generalized 8-vertex 
model. The next step is to anti-symmetrize the action by subtracting from 
the above expression (0) the same terms but with reversed order of the 
Grassmann variables and dividing by 2. It is convenient to denote the re- 
sulting action using matrix notation. We order the Grassmann variables in 
a vector (T denotes transposition) 



rj(x) = ^ + i(x),r/_i(x),r/ +2 (x),r?_ 2 (a;)) 



and define the 4x4 matrices P±i and P±2 



(6) 



P+l(i,j) 



They obey P l 



■^2^,3^,4 



P-l(i,j) 



-P-y. We also define 
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We remark that for writing in this form the free fermion condition (||) 
was used. The determinant of \i is given by det \i = w\ and pL T = —fi. 
Using the matrix notation the action reads 

%] = \^2v{x) T K{x,y)r](y) , 
with the kernel K consisting of 2 parts K = —M + R, where 

±2 

M(x,y) = jj, 5 XtV , R{x,y) = ^ B v {x)P v S x+i>j y , 

u=±l 

and we have defined 



B-„(x) = B u (x-v) , v = 1,2. 

It has to be stressed, that here the link variable is not complex conjugated 
when hopping in negative ^-direction. This is different from lattice gauge 
theory (see e.g. |T^]) where the link variables are conjugated for backward 
hopping. This difference indicates that for lattice field theories with fermions 
coupled to a vector field a different structure emerges. 

The partition function can now be written as a Pfaffian, and since by 
construction the kernel K of the action is anti-symmetric, the Pfaffian is 
given by the root of the determinant of K. An overall factor det[— M] = w\ 
can be extracted and the rest is expanded in powers of the hopping matrix 
giving 

Z = (-1)I A I J Dr) a*"**" = (-1)I A I Pf K 

= (-1) |A| ^/det[-M + R] 

= (-l) |A| ^det[-M]Y/det[l - M~ l R] 

= (-Ijl^e-^^, 

(8) 
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where we introduced the hopping matrix 

1 ±2 

H(x,y) = M- l R{x,y) = — V B„{x)T v <W, y , (9) 
with IV = wi^i^Py, or explicitly 

\ 


/ 

+w 6 \ 
-w$ 
+w 3 
J 

The series in the exponent of (@) converges for ||-£f|| < 1. By choosing 
suitable values for the weights Wi and the fields B fM (x) this can be achieved. 
Note however, that this restriction is only a technical assumption due to the 
expansion used here. The final result is simply a finite polynomial in the Wi 
and Bn(x) and the restriction of the variables can be lifted then. 

Due to the Kronecker deltas in (|9|) the traces TrH n reduce to products 
of matrices T u supported on closed loops I. Since on the rectangular lattice 
all closed loops have even length, traces over odd powers of H vanish. For 
even powers we obtain 




By Cr> we denote the set of all closed, connected loops of length 2k and 
base point x. The first product picks up a factor B v (y) whenever the loop 
I G C^ k ^ runs through the link (y,f). The last term is simply the trace 
of the product of matrices as they appear along the loop. By explicit 
calculation one can show that 
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Y±n T Tfl = , n = 1,2 . 

This implies that whenever the loop I turns around at a site and runs back 
on its last link the contribution of this back-tracking loop vanishes. Thus the 
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set C x only contains closed, connected, non back-tracking loops of length 
2k with base point x. For each loop both possible orientations have to be 
taken into account. The loops in C x 2 ^ can self-intersect and retrace parts 
of their contour or even their whole contour. 



4 Computing the traces and final result 

In this section we evaluate the traces over the matrices T^. The loops 
are decomposed into 4 simple elements (depicted in Fig. and 
algebraic relations for products of corresponding to these elements can be 
used to give a constructive procedure for determining the traces TrFJ^ei 
(compare also [20] where the case of the Ising model in this expansion was 



discussed in more detail). 

The first of the four elements, basic loop, is the closed loop around a 
single plaquette (compare Fig. ||.a). One finds by direct evaluation 

TV r^r^r^ = -w 5 w e w 7 w 8 , (11) 

where /i, v = ±1,±2 , n ^ ±za Thus for the basic loop around a single 
plaquette we always obtain —1 times the product of the weights for the 4 
tiles showing a corner (compare Fig. |]), independent of the starting point 
and the orientation. 

The second element is referred to as telescope rule. The telescope rule 
can also be shown by explicit evaluation and reads 

r±ir ± i = w 4 r ± i , r ±2 r ±2 = w 3 r± 2 . 

The geometrical interpretation of these two identities is simple. A piece 
of loop can be shrunk or stretched in horizontal (vertical) direction and a 
factor tt>4 (W3) has to be collected (compare the graphical representation in 
Fig. |.b). 

The third rule is the kink rule depicted in Fig. ^.c. Algebraically it reads 

r ±M r ±I ,r ±At = w 5 wq r± u , n, v = 1, 2 , n ^ v , 
Fi^r^r-^ = w 7 w 8 r± u , n, v = 1, 2 , ^ ^ v . 

This rule, together with the telescope rule allows to remove superfluous kinks 
or corners in a loop (or sub-loop) and to reduce this loop (sub-loop) to rect- 
angular form which then can be shrunk to a loop (sub-loop) around a single 
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= - w w w w 



= - w w w w 

3 4 7 8 



+ rotated graphs 



Figure 2: The four elements used for the computation of the traces: a: 
basic loop, b: telescope rule, c: kink rule, d: intersection rule. 



plaquette using the telescope rule. One simply replaces a piece of loop by 
a simpler piece (i.e. one with less corners) which can be constructed using 
telescope and kink rule (see [^0| for a detailed description of the procedure) . 
During this process one has to collect all emerging factors Wi. In case the 
loop we started with contained no sub-loop, this procedure will reduce it to 



the basic loop and using (11) the result for the trace can be read off. 



In case the loop contains sub-loops, we first reduce a sub-loop which has 
only one self-intersection (such a loop always exists) to a sub-loop around 
a single plaquette. The sub-loop can then be removed by the fourth rule, 
the intersection rule. It states, that when removing a sub-loop, in addition 
to the corner and monomer factors, a factor of —1 has to be collected. The 
algebraic expression reads (for the graphical representation see Fig. |2|.d) 



-w 3 w4w 5 w 6 r ±M r Ti , , 



for fj,, v = 1, 2 , n ^ v. As outlined above, and discussed in more detail in 
TO], the four elements can be combined to compute the trace in a construc- 
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tive way. After collecting all monomer, corner and intersection factors, we 
find the result 

^ n r M = -(-i) sw n< i( °. (12) 

/.tGi i=3 

Here s(l) denotes the number of self- intersections the loop / has. The ex- 
ponents rii(l) give the numbers for the abundance of the line elements as 
they are depicted in Fig. |l]. E.g. when the loop changes from heading east 
to heading north at a site, it picks up a factor of w$, and similarly for the 
other tiles. Note that the loops I in the hopping expansion appear as an 
ordered set of instructions for the directions the loop takes as it hops from 
one site to the next. The element corresponding to tile Nr. 2 with weight 
u>2 (compare Fig. [I]) does not appear. We also remark, that the result ( |T2| ) 
is independent of the orientation of the loop. 
Inserting (12) and ( ^0| ) into (|8|) one finds 




Z = (-i)iM A| expf^E^E e (-i) s(0 n My) n- m(/) 



(13) 

The final step is to eliminate the explicit summation over the base points. 
When doing so a minor subtlety has to be discussed. Consider first a loop 
which runs through its contour only once. Each of the 2k lattice points it 
visits can serve as the base point producing a factor 2k which cancels the 
corresponding factor in (|l3|). If now a loop I completely iterates its contour 
/(Z)-times, then there are only 2k /I (I) different choices for a base point and 
a factor 1/1(1) remains. The final result is 

Z = (-l)'M A| xexp(£i^(J-f n B v (x) f[w?A. (14) 

\iec W Wl (x,v)ei i=3 J 

Here the sum runs over the set C of all closed, non back-tracking loops of 
arbitrary length. Each loop is included with only one of its two possible 
orientations leading to the cancellation of the factor 1 /2 which appeared in 
the last expression. By |/| we denote the length of the loop I and 1(1) is the 
number of iterations of its complete contour. We remark, that (|i~4] ) does not 
explicitly contain W2 , but this weight is related to the other weights through 
the free fermion condition (M). 
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5 Discussion of the result 



Several aspects of the result ( |i~4| ) should be discussed. When viewing vertex 
models as models of Peierls contours for Ising-type spin models the exponent 
in (|l4]) is proportional to the free energy in terms of loops. For the standard 
Ising model a representation of this type is given in [21]. Equation ( |T4| ) 



generalizes this results to 8- vertex models coupled to an external field which 
also includes the case of Z^-spi 11 models with locally varying couplings. 



The expression (14) is also remarkable from an algebraic point of view. 



In principle, it is possible to expand the exponent in (14) in a power series of 
w\ ~ 2k . The coefficients in this series can be expressed in terms of products of 
loops / £ C. In turn one would like to identify these loops with the contours 
C(t) which occur for the tilings t in the original formulation of the model. 
The loops I £ C, however, can occupy links several times and also taking 
products of these loops will create multiply occupied links while the contours 
C(t) only have empty or singly occupied links. When analyzing some terms 
of this expansion, one finds however, that due to the self-intersection factor 
(— 1) S W all contributions where links are multiply occupied get cancelled. 
Only terms without multiply occupied links survive and these correspond to 



the contors C(t). Our equation (14) establishes a remarkable feature of this 
expansion, namely that the cancellation of terms with multiply occupied 
links is independent of the corner and monomer weights. When applying 
the formula to fermionic lattice field theories (see below) this cancellation of 
all terms with multiply occupied links implies a considerable simplification 
of the hopping expansion for the fermion determinant. 

Probably the most interesting application of Eq. ( |l4| ) is its use when 
relating fermionic lattice field theories to vertex models. An example of 
this type of application of our formula Eq. (|l4|) is discussed in Q, where 
it is proven, that 2-D Wilson fermions in an external scalar field <f(x) are 
equivalent to our model in its form (|2|) with w\ = l,u>2 = 0,^3 = W4 = 
k,W4 = ... u>8 = n/y/2, where k is related to the mass parameter through 
k = (2 + m) -1 . For the case of free Wilson fermions {<p{x) = 1) this was first 
shown by Scharnhorst [p3] with different methods (for a presentation using 



the techniques of this paper see [24]). [|| extends this result to non-trivial 
background fields using a careful analysis of the hopping expansion for the 
Wilson fermions. Wilson fermions give rise to a bilinear Grassmann action 
(instead of a quadratic form here) and the object to be expanded is the 
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fermion determinant (instead of the Pfaffian). The Dirac 7-matrices play 
the role of the hopping generators which we encounter here, and also for 
the Wilson fermions it is possible to compute all the traces which appear 
in the hopping expansion. The final step in the proof of equivalence is the 
identification of the correct weights u>i by comparing the hopping expansions 
for the two models. The details of this mapping are given in 0. Further re- 
sults are obtained in Q by integrating out the scalar field, which can e.g. be 
used to generate the Gross-Neveu model and map it to another vertex model 
(for similar mappings of the Schwinger model and the 2-D Thirring model 
with Wilson fermions see [^, [25|, |26| ) . Such mappings of lattice field theo- 
ries onto vertex models (loop representation) are powerful tools, since often 
the numerical simulation of the vertex model is simpler than analyzing the 
original model and allows for higher precision (see e.g. [^J where this was 
demonstrated for the case of the strongly coupled Schwinger model). In par- 
ticular the loop expansion of the fermion determinant developed here and 
in Q overcomes the fermion sign problem. This permutation sign plagues 
fermionic systems in the Hamiltonian approach |^8| as well as in the Grass- 
mann path integral formulation |2{| and in a numerical treatment requires 
an exponentially increasing number of Monte Carlo configurations as the 



volume or inverse temperature are increased [28]. Currently we are working 



on an implementation of new numerical algorithms which are based directly 



on the loop representation |30fl . 
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